Effect of interactions on the cellular uptake of nanoparticles 
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We present a simple two-state model to understand the size-dependent endocytosis of nanopar- 
ticles. Using this model, we elucidate the relevant energy terms required to understand the size- 
dependent uptake mechanism and verify it by correctly predicting the behavior at large and small 
particle sizes. In the absence of interactions between the nanoparticles we observe an asymmetric 
distribution of sizes with maximum uptake at intermediate sizes and a minimum size cut-off below 
which there can be no endocytosis. Including the effect of interactions in our model has remark- 
able effects on the uptake characteristics. Attractive interactions shift the minimum size cut-off 
and increase the optimal uptake while repulsive interactions make the distribution more symmetric 
lowering the optimal uptake. 



INTRODUCTION 

The endocytic process [l], is of paramount importance 
to understanding the cellular uptake of nano-materials, 
essential for the development of gene and targeted drug 
delivery tools. A key feature in the development of such 
tools is to achieve effective cytosolic delivery. To this 
purpose there have been experiments using liposomes Q , 
nanoparticles (NPs) [4 [jpolymerosomes [1431, nanotubes 
d, 0] , electroporation l3| and ultrasonic treatments [ll| . 
Several of these techniques also suffer from the problem 
of high levels of cytotoxicity although recent experiments 
using polymerosomes overcome this shortcoming. 

In several of these experiments using gold and silver 
nanoparticles, nanotubes and polymerosomes [sI-ItI, [ol [l2|- 
17[ particle size plays an important role in the cellular 
uptake. These experiments suggest that endocytosis of 
NPs is receptor-mediated and that there is an optimal 
size where the uptake is maximum. Most theoretical ap- 
proaches IB" 21 1 to study the effect of NP geometry on 
cellular uptake predict a threshold radius below which 
there can be no cellular uptake, and an asymmetric distri- 
bution of the uptake which decays with particle size. Al- 
though these approaches correctly predict the size where 
the uptake is optimal (~ 20 — 30 nm), experimentally the 
distribution is symmetric and does not seem to agree with 
the lower bound as predicted from the theories [l^ [3 ■ 
The answer to this anomaly could be hidden in the highly 
complex endocytic mechanism itself. 

The endocytic process involves the selection and segre- 
gation of the cargo at the cell surface, subsequent invagi- 
nation and pinching off from the cell membrane, and, fi- 
nally, the transport of these vesicles to intracellular com- 
partments where they fuse with the target membrane. 
The mechanisms by which specific cargo are internal- 
ized differ in their morphological and biochemical details 
P, 0, |2l|- However, recent evidence suggests the 
need to consider the sharing of molecular machinery de- 



pending on the nature of the cargo and to understand 
the basic physical principles common to these different 
uptake mechanisms. 

A key step in the endocytic process is the segrega- 
tion and clustering of cargo on the cell membrane which 
are believed to be the sites where molecular machin- 
ery may be recruited to generate membrane curvature, 
form membrane invaginations and subsequently cause 
scission P, [2^ . The mechanisms of formation of these 
nanodomains can be both passive and active. Although 
passive clustering which does not involve ATP hydrolysis 
can occur via intermolecular cargo interactions, Reynwar 
et. al. ^231 showed using coarse grained simulations, that 
curvature-inducing model proteins adsorbed on lipid bi- 
layer membranes could experience attractive interactions 
that occur purely as a result of membrane curvature. 
These interactions could result in clustering and subse- 
quent invaginations. Thus passive clustering can happen 
even in the absence of specific cargo interactions. An 
example of the passive clustering is the binding of Shiga 
toxin — a bacterial toxin — to glycolipid receptors, Gb3, in 
the cell membrane of certain cell types [2^, [2^. Although 
the Shiga toxin molecules do not interact directly, they 
induce the clustering of the Gb3 lipids, thereby causing 
local membrane curvature. 

The active mechanisms for cell surface clustering re- 
quires ATP and are therefore energy dependent pro- 
cesses. An example of active clustering is the formation 
of nanoscale clusters of glycosylphosphatidylinositol an- 
chored proteins (GPI-AP) on the cell membrane which 
are required for their subsequent endocytosis [27| . These 
proteins exist as monomers and nanoclusters ('^ 4 — 6 nm 
in size, consisting of < 5 molecules) with the interconver- 
sion between the two being spatially heterogeneous, being 
coupled to an active cortical cytoskeleton. Perturbing the 
cortical actin activity affects the construction, dynam- 
ics and spatial organization of these nanoclusters [28| . 
This type of active segregation also happens with gan- 
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glisides GMl and GM3 in the exoplasmic (outer) leaflet 
|29l] and Ras isoforms in the cytoplasmic (inner) leaflet 
of the plasma membrane (soj . 

The above examples of passive and active clustering 
do not involve the clathrin mediated endocytic pathway 
where the entrapment of the cargo occurs by its associa- 
tion with adapter proteins. However, even in the clathrin 
coat mediated endocytic pathway, the clathrin lattice 
organizes the epsins or BAR domain proteins into do- 



mains, which then locally deform the membrane [311.13 



Therefore the segregation and clustering of cargo on the 
cell surface is highly important in the endocytic process. 
This naturally raises the following questions: could cell 
surface clustering affect the size-dependent cellular up- 
take of NPs and if so how could we model the clustering 
process? Recent experiments have shown evidence 
of NP surface clustering on the cell membrane and it is 
important to investigate this in some detail. 

Here, we study systematically for the flrst time, the 
effect of interactions on the cellular uptake of NPs us- 
ing a thermodynamic model first proposed by Tzlil et al. 
[3, 111] and subsequently studied by Zhang et al. [l^l in 
the NP context. We develop our model by incorporating 
interactions between NPs. Using a simplified two-state 
version of the model, we then elucidate the relevant en- 
ergy terms which affect the uptake of NPs in the absence 
of interactions. We then show that interactions between 
NPs indeed affect the minimum radius of uptake as well 
as the distribution. 



THE MODEL 

In this model, the system, which consists of a cell and 
ligand-coated spherical NPs in a solution, is in a thermo- 
dynamic equilibrium at which a certain number of NPs 
are endocytosed. At this state, N NPs adhere to the 
cell surface that contains L receptors via ligand-receptor 
binding, and are wrapped to different extents by the cell 
membrane (see Fig. [1}. The receptors diffuse freely on 
the cell surface and are segregated into Lp free recep- 
tors in the planar membrane and Lf, bound receptors in 
the curved regions. Let A be the cross-sectional area 
of a receptor. For a given NP radius R, the number of 
receptors that can attach to the NP is K — AttR'^/A. 
For convenience, we shall choose A as our unit of area 
and VA = R^Att/K as our unit of length. The total 
membrane area is denoted by MA (in units of A) and is 
therefore the total number of sites on the membrane that 
are accessible to the receptors. The surface concentration 
of NPs is then c = N/MA. 

We have assumed that the time scale for endocytosis is 
much larger than the time for the receptors to diffuse and 
segregate into the curved and planar regions. Therefore 
the distribution of wrapping sizes and receptor densities 
can be treated using equilibrium statistical mechanics. 



Let Uk denote the number of NPs wrapped by a mem- 
brane section of area k where k varies discretely be- 
tween fc = (unwrapped state) and k — K (completely 
wrapped state). Then, we have 
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k=0 



where Mi,A is the total curved membrane area associ- 
ated with the wrapped NPs and MpA = (M — Mi,)A is 
the total area of the planar regions. The binding of a 
ligand and a receptor releases chemical energy, e, which 
drives the wrapping at the cost of the energy required to 
bend the membrane. Therefore, the diffusion of recep- 
tors inside the curved regions should lower the energy of 
the system facilitating wrapping. However, this leads to 
the segregation of receptors between planar and curved 
regions, which costs entropy. Also diffusion of free re- 
ceptors into the curved regions increases the total curved 
area (more ligand-receptor bonds) and hence increases 
the total membrane bending energy. Furthermore, at- 
tractive (repulsive) interactions between the NPs could 
lead to clustering (anti-clustering) and therefore affect 
the wrapping size distribution of NPs. To determine the 
size distribution of the varyingly wrapped NPs, we first 
write down the free energy of the system as. 
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= Mp [(j)p In 0p + (1 - (j)p) ln(l - <j)p)] 

+ Mb [<pb In 0b + (1 - (f>b) ln(l - 0b)] 
+ nk[ln (nk/M) - 1] - eU + kMb 

k 

+ ^nfeAfc +^nfcrfc + w^fcfc'nfenfe/, (3) 



where (pp = Lp/MpA and <j)b = Lb/MbA = {L-Lp)/MbA 
denote the densities of the receptors in the planar mem- 
brane and the wrapped regions respectively, ks denoting 
the Boltzmann constant and T being the temperature. 

The first three terms in the free energy are entropic 
contributions written in terms of a two dimensional lat- 
tice gas model: 

• Alp J0p hi(j>p + {l — (pp) ln(l — (j)p)] represents the con- 
figurational entropy of Lp free receptors distributed 
in the Mp sites of the planar parts of the membrane. 

• Mb[0b In 0b + (1 — 0b)ln(l — 0b)] represents the 
configurational entropy of distributing Lb receptors 
among the Mb sites of the curved regions. 




FIG. 1: Schematic figure of adliering NPs wrapped by cell membrane. The NPs are wrapped to different degrees with some 
of them being internalized, (a) Single NP wrapping showing ligand-receptor binding, (b) Cluster of NPs wrapped by the cell 
membrane. 



• ^ J, rife [In (rife /M) — 1] is the configurational entropy 
of a 2D mixture of wrapped NPs when treated as a 
multicomponent ideal gas. However, we are inter- 
ested in interacting NPs and would therefore have 
to include interaction energy for this 2D mixture. 

The next five terms are energetic: 

• ^eif, — — A/f,0be — ~(t>b^J2k ^"•fc ^^^'^ total chem- 
ical energy released upon the binding of Li, ligand- 
receptor pairs. 

• kMb = Kj2k ^'^fc is the total membrane curvature 
energy in the budding regions. For a spherical ge- 
ometry, the bending energy per unit area across 
a NP of radius R is k ^ {KA/2kBT){2/ R - co)^ 
where k denotes the bending modulus. Note that 
the spontaneous curvature (cq) of cell membranes 
is nonzero. In our analysis we consider a van- 
ishing spontaneous curvature (cq = 0). There- 
fore, k = 2KA/kBTR'^ = STTn/kBTK. A positive 
spontaneous curvature (cq > 0) could decrease the 
value of k while a negative spontaneous curvature 
(co > 0) could increase the value of k, for a given 
value of K. We shall discuss the effect of sponta- 
neous curvature on cellular uptake later. 

• ^i^nkTk'. Total work of pulling excess membrane 
towards the wrapping sites against lateral tension 
a. For a single NP wrapping, the excess area 
pulled towards the wrapping site is AnR^k^ / — 
k'^A/K [s^]. Thus the excess energy is F^ = 
a X excess area — k^aA/kBTK . 

• nfcAfc : Total line energy of the rim, where A(fc) 
denotes the line energy of a fc-bud. Assuming a 
spherical shape of the membrane at the rim of a 



partially wrapped NP, A^ is modeled as being pro- 
portional to the length, Ck , of its rim, with a con- 
stant line energy per unit length 7 [18|. Therefore, 



(4) 



Note that Ck vanishes for A; = and k — K and is 
maximum {2ttR) for a half- wrapped NP (fc — K/2). 
However, the local wrapping behavior of the mem- 
brane to a NP is different from the assumption 
made above (s^, 35 1. We need to consider an ad- 
ditional bending energy for the unadsorbed mem- 
brane detaching from the NP at the rim. Although 
the fc— dependence of this energy is different from 
the simple form assumed in Eq. (|4]), the general 
features of large energies for half-wrapped state 
(fc — K/2) and very small energies for unwrapped 
and completely wrapped states are the same. 
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is the interaction energy between the partially 
wrapped NPs, w denoting the strength of the inter- 
action (or the second virial coefficient). We assume 
the interaction to depend on the degree of defor- 
mation and curvature of nearby curved membrane 
patches and therefore to the degree of wrapping of 
the cell membrane to individual NPs. Thus the to- 
tal interaction energy when summed over is propor- 
tional to the total curved area. Interaction between 
membrane inclusions or adsorbates could arise due 
to a variety of different mechanisms (36l | including 



membrane fluctuations 
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In the final stages of endocytosis the membrane wrapped 
NP pinches off which results in a topology change. This 
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FIG. 2: Size dependent cellular uptake of non-interacting NPs. (a) The red line is for the two-state model with = = r^, 
the green (□) and blue (A) points are for the _ftr-state system with Ffc = and Tk respectively, (b) The red line is for the 
two-state model while the blue (•) points indicate the large R behaviour. 



severing mechanism of the wrapped NP from the mem- 
brane is brought about by proteins such as dynamin 
and C-terminal binding protein 3/brefeldin A-ribosylated 
substrate (CtBP3/BARS). According to Gauss-Bonnet 
theorem '39[, this leads to an increase of inn in Gaus- 
sian bending energy, with R representing the Gaussian 
bending rigidity of the cell membrane. In our model we 
are interested in events prior to the final scission process 
and will therefore safely ignore this term. In fact, all par- 
tially wrapped NPs for which k > 0.9K will be assumed 
to be endocytosed. 

To find the equilibrium state of the system we mini- 
mize the free energy with respect to Li, and rik- From 
dF /dL\j = 0, we have 



(5) 
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Minimizing F subject to the constraint of Eq. ([1]), we 
get the normalized wrapping size distribution as 



(6) 



where we have defined 

a = e^-«, (7) 

and 

Pk = Afe -t- Tk + '2.wk ^ k'nik' 

k' 

= Ak+Tk+2wcMkJ2k'Pk'- (8) 

k' 

The conservation condition for the receptors gives 

4ip{l - c^kpk) + (pbc^kpk ^ (t)o. (9) 



(/>p and back into Eq. ^ yields the wrapping size 
distribution, and hence the number of fully internalized 
NPs 



riK = cMpK- 



(10) 



Therefore, uk gives the cellular uptake of nanoparticles. 
We first study the effect of particle size on the cellular 
uptake of non-interacting NPs and then see the effect of 
interactions on the distribution. 

The choice of the physical constants is mostly guided 
by experimental data although there are some free pa- 
rameters as well. The bending modulus (k) of biomem- 



branes is typically on the order of 20 k^T [40|, |41 [. 
The receptor-ligand binding energy, e is assumed to be 
comparable to antibody-antigen interaction and is esti- 
mated to be on the order of 15 - 25 fceT [H, The 
length scale is set by the length of the receptor (VA) 
which is typically on the order of 15 nm. Therefore, 
A ~ 225 nm^. Experimental information suggests that 
the number of receptors varies from 50 — 500 per /im^ 
E El, El El- This implies that 00 could vary from 
0.01 to 0.1. The concentration of NPs, c, can vary be- 
tween 0.001 and 0.005. The diameter of the cell being 
« 15 nm, the surface area of the cell is « 707/im^. There- 
fore, M = 3.14 X 10^. In our numerical analysis, we 
choose K = 20 feB^, e = 25 /cbT, c = 0.003, (/)o = 0.05, 
and M = 3.14 x 10^ [Ulli]. In what follows, we choose 
a — and consider the effect of a on uptake in a later 
section. Both 7 and w are free variables and we choose 
7 = 1.0 (in units of ksT per unit length, ^/A). w is 
varied from zero (non-interacting) to positive (repulsion) 
and negative (attraction) values. 



TWO-STATE MODEL 



The densities of the receptors 4>p and 4>i, can be obtained 
by numerically solving Equations ([S])-®. Substituting 



To analyze the size-dependent uptake of NPs we make a 
major simplification in the model. We assume that the 
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NPs upon arrival to the cell surface are either endocy- 
tosed completely or remain free without there being any 
intermediate wrapped state. Then the model essentially 
reduces to a two-state model with the two states being 
k = and k = K. Our goal is to come up with the 
minimal model to understand the experimental uptake 
behavior and to find the only relevant energetic contribu- 
tions. Note that the line energy term (Afc) automatically 
vanishes with this simplifying assumption. We now study 
the uptake behavior (i) in the absence of interactions and 
(ii) when the NPs interact. 



(i) Non-interacting case (w — 0) 

With the two-state model, in the absence of interactions 
and with A^^ = = Ffe, we observe [Fig. [^fa)] that be- 
low a critical radius, i?minj there is hardly any uptake. 
Above i?min7 the uptake increases sharply to reach a max- 
imum and then decays as a power law with increasing 
radius. Thus, the two-state model correctly reproduces 
the optimal uptake behavior at intermediate radii seen in 
experiments. We compare our results for the two-state 
model with the full K-state model both in the presence 
and absence of the Tk term. We find that for 7 = 1 and 
a = 0.001 (in units of fcsT per unit area) the surface 
tension term does not affect the size-dependent distribu- 
tion significantly. As we shall show later, increasing a 
decreases uptake significantly although the uptake char- 
acteristics remain unchanged. Thus we conclude that 
the only relevant energy terms in understanding the size- 
dependent endocytosis of NPs are the energy released on 
ligand-receptor binding and the energy cost in bending 
the membrane. We verify this further by analyzing the 
behavior at large and small radii. 

Behavior at large R. For the two state model, Eq. ^ 
(the conservation condition) reduces to giving the com- 
pletely wrapped particle size distribution as 



Pk 



1 



(11) 



In the large R limit (or large K limit), the receptor 
density in the fully enveloped NPs is almost saturated, 
(ph ^ 1, whereas the free receptor density is negligible 
6„ ~ 0. Then 



Pk 



cK 



lA 



Therefore, cellular uptake for larger NPs is inversely pro- 
portional to the square of the size of the NPs and repro- 
duces the numerically predicted behavior at large radius 
exactly [Fig. EJb)]. 

Behavior at small R. To understand the low cellular 
uptake at smaller radii we follow Tzlil et al. to make 



the macroscopic (bud) phase approximation, i.e. assume 
that instead of the curved regions being made up of sev- 
eral NPs wrapped to difi^erent extents, there is a single 
NP with wrapped area Mf, that coexists with the planar 
membrane phase. This approximation causes the config- 
urational entropy of the NPs in the free energy expression 
to drop off. Also A^ = for all k. Minimizing the re- 
sulting free energy with respect to Lb gives Eq. ([5|) and 
minimizing it with respect to Mb yields 



1 



1 



1 



l~0fc 



(13) 



Solving Eqs. ([5]) and p3)) we can determine the receptor 
densities in the two coexisting phases as [3] 



1 



and (f>p 



1 



(14) 



Therefore, we can have coexistence between the planar 
and wrapped phases only if e > k > 0. Thus, for a single 
wrapped NP, e = k is the critical value below which we 
cannot have wrapping. Substituting for k, we get the 
critical radius for the onset of wrapping as [46|, 



^y2^^A/e 



(15) 



For the values of k and e used in the numerical estimates 
we get i?min ~ 19 nm. 

(ii) Interacting case {w ^ 0) 

As observed above, the uptake of NPs in the absence 
of interactions is highly asymmetric and also predicts 
a lower radius cut-off. Experimentally, the distribution 
has found to be rather symmetric both for Au nanopar- 
ticles [H, [1^ and DNA wrapped single-walled carbon 
nanotubes (DNA-SWNT) 0. Moreover, there is sig- 
nificant internalization of particles below the minimum 
radius predicted by the model. Experiments using DNA- 
SWNT show an increase in near-infrared fluorescence 
from SWNT concentrated at the external cell membrane 
during the early stages of endocytosis mechanism Q , in- 
dicating possible clustering of nanotubes on the cell sur- 
face prior to uptake. To incorporate clustering in our 
model, we include an effective interaction in our model. 
The idea is that interactions could lead to clustering 
which could drive wrapping of NPs of smaller sizes. In 
our model, the interaction between NPs is controlled by 
the interaction parameter, w. Negative w implies at- 
traction while positive w implies repulsion. We do our 
analysis for the two-state model. 

Results for attraction {w < 0). Figure [3Ka) shows 
the results for size-dependent uptake of NPs in the pres- 
ence of attractive interaction. Similar to the behavior 
in the absence of interactions, we flnd that below i?min, 
there is hardly any endocytosis. Above Rmim the uptake 
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FIG. 3: Size dependent cellular uptake of interacting NPs. (a) The curves are for different values oi w = 0.0 (red, □), 
w = —0.00002 (green, o), w = —0.00003 (blue. A), and w = —0.00004 (pink, o). -Rmin decreases with increasing strength while 
the optimal uptake increases, (b) The curves are for different values of u; = 0.0 (red, □),«; = 0.00001 (green, o), w = 0.00005 
(blue. A), and w = 0.0001 (pink, o). Optimal uptake decreases with increasing repulsion. 



increases rapidly, and subsequently reaches a maximum 
and then decays slowly. We note, however, that the up- 
take mechanism is strongly dependent on the value of w. 
With the increase in the strength of the attractive inter- 
action, (increasing the minimum radius for uptake, 
-Rmin, decreases substantially. Also the maximum uptake 
increases with increasing \w\ indicating that uptake be- 
comes more favorable in such circumstances. To explain 
this behavior we need to look carefully at the energet- 
ics. Without interactions, the low radius cut-off is deter- 
mined by the competition between the energy released on 
ligand-receptor binding and the energy cost in bending 
the membrane. Attractive interactions cluster NPs and 
lower the free energy locally so that the system can afford 
to pay the cost of membrane bending energy, resulting in 
endocytosis for smaller particle sizes that would other- 
wise have been prohibited. 

Results for repulsion {w > 0). Figure [3l^b) shows 
the results for size-dependent uptake of NPs in the pres- 
ence of repulsive interaction. We find that the behavior 
is distinctly different from that with attractive interac- 
tions. Although the uptake mechanism strongly depends 
on the interaction strength w, the characteristics differ 
significantly. Interestingly, the lower radius cut-off i^min 
does not shift with increasing w and is the same for all 
w values. The maximum uptake however decreases with 
increasing w. The uptake behavior at large R is also af- 
fected. Although the uptake decreases at large R, the 
decay is much slower as w increases. This behavior can 
again be explained by looking at the energetics. Repul- 
sion between NPs (positive w), pushes them apart and 
the energetics at lower radius are governed by the same 
single NP wrapping energetics as in the absence of inter- 
actions. Therefore, i?min does not change when we in- 



crease w. However, the uptake for R > i?inin is affected 
since positive w increases the global free energy cost, thus 
decreasing the optimal uptake and affecting the behavior 
at large radius making the uptake more symmetric with 
particle size. 

Therefore, we find that interactions affect cellular up- 
take of NPs significantly. We studied the two cases of 
attractive and repulsive interactions separately. How- 
ever, in the physical system, we expect both of these 
interactions to be present simultaneously with attractive 
interaction affecting the uptake of smaller particles and 
repulsion dominating for the bigger ones. Such a picture 
could provide a qualitative understanding of the exper- 
imental observation of size-dependent cellular uptake of 
NPs. 



Effect of membrane tension and spontaneous 
curvature 

The effect of membrane tension is to lower the optimal 
uptake of NPs 47, 4^. In Fig. |4]we show the variation of 
uptake with system size for different values of the surface 
tension a (in units of ksT per unit area) in the absence 
of interactions. The nature of cellular uptake with the 
particle size remains the same for different values of a 
although the amount of uptake decreases. 

In all our analysis we have ignored the effect of spon- 
taneous curvature (cq). As we have mentioned earlier, 
a non-zero cq could lower or raise the value of k which 
should in turn lower or raise the free energy barrier for 
the uptake of a nanoparticle. Hence the optimal uptake 
could indeed be lower or higher depending on the value 
of Cq. However, the characteristics of uptake behavior is 
not altered. 



7 



6000 




Nanoparticle Radius (in nm) 

FIG. 4: Size dependent cellular uptake of non-interacting NPs 
for different values of a = 0.001 (red, □), a = 0.02 (green, o) 
and CT = 0.1 (blue. A). 

SUMMARY AND CONCLUSION 

In this paper, we have studied the effect of interactions 
between surface-bound NPs on their subsequent endocy- 
tosis in the context of a statistical thermodynamic model. 
One of our first results is to clearly show the relevant 
energy terms required to understand the uptake mecha- 
nism qualitatively. To do so, we simplified our model to 
a two-state model, which captures the essential features 
of the uptake behavior. We showed that apart from the 
important entropic contributions coming from the distri- 
bution of the receptors and the NPs on the cell surface, 
the energy terms that dictate the uptake characteristics 
are the energy released via ligand-receptor binding and 
the energy cost in bending the membrane. Although we 
show that the line energy and energy in pulling excess 
membrane area are not significant when the aim is to un- 
derstand the specific uptake behavior, we do not rule out 
the importance of these terms in the endocytic process. 
These energy terms, as well as the Gaussian bending en- 
ergy during the pinch-off, are relevant for the studies of 
membrane invaginations and wrapping. However, they 
may be dispensable when addressing the question of the 
number of NPs endocytosed at a given time. 

Beyond the two-state model, our results show that in- 
teractions between NPs could have a drastic effect on the 
uptake process. Attractive interactions lead to clustering 
of NPs, which effectively lowers the free energy threshold 
for wrapping and therefore shifts the lower cut-off radius. 
This is not possible in the absence of interactions un- 
less, of course, we changed the relative values of k and e. 
For fixed k and e, we see that repulsive interactions also 
have a significant impact, in that they cause the cellular 
uptake to be reduced and modify uptake characteristics 
towards being more symmetric. 

In our model, interactions are assumed to be propor- 
tional to the wrapped area and the strength of the in- 



teraction, w, is varied freely. Therefore, w could be 
thought of as an effective parameter which models the 
membrane- mediated forces [s^. Miiller et al. [i^ [s^l 
study the curvature-mediated interaction between parti- 
cles on the cell surface. From a purely geometric analysis, 
they showed that the net force between the particles is 
due to a competition between the force associated with 
the curvature along the direction joining the particles 
(which leads to repulsion) and the force associated with 
the curvatures perpendicular to it (which leads to attrac- 
tion). It will be interesting to use this kind of model for 
membrane-mediated interaction in a more detailed study 
of receptor-mediated endocytosis of interacting NPs. 
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